Reply to the Comments on the XENONIOO first Dark Matter Results 

The XENONIOO Collaboration 
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The recently submitted preprint on the first results 
from the XENONIOO dark matter experiment [1] was fol- 
lowed by a criticism by J.I. Collar and D.N. McKinsey [2], 
focused on our extrapolation of the scintillation efficiency 
£eff to the lowest nuclear recoil energies, where no data 
and no theoretical model exist. Here we add clarifications 
on our analysis and comment on their criticism. 

The main XENONIOO result demonstrates, with only 
11 days of data, the potential of this ultra-low back- 
ground experiment to exclude new parameter space in 
spin-independent WIMP-nucleon cross section. The high 
mass WIMP limit is unaffected by the controversy on 
the value of i2eff- The light WIMP interpretation of the 
DAMA [3] and CoGeNT [4] data is excluded at 90% CL 
for a constant extrapolation of the global fit of all >Ceff 
data from fixed-energy neutron experiments (see Fig. 1 of 
our preprint). However, when >Ceff is assumed to follow 
the 90% lower limit of the global fit with a conservative 
logarithmic extrapolation to zero near 1 keV^, a fraction 
of the CoGeNT parameter space remains uncovered. 

We further argue that we don't see any events even 
down to a threshold of 3 PE, departing from our a priori 
chosen threshold of 4 PE, where, however, our acceptance 
is still quite high (as shown in Fig. 2 of the preprint). 
This excludes all of the DAMA favored region, with ion 
channeling, at 90% CL, even for the conservative case 
of £eff- In this case, we would expect one event for a 7 
GeV/c^ WIMP, at the lower edge of the CoGeNT favored 
region. With our current data, only this conservative 
Cef[ leaves room for low mass WIMPs compatible with 
CoGeNT in the mass range of ~7-9 GeV/c^. 

We will now comment on some of the statements made 
in [2] in more detail. 

1. In [2] it is stated that we use a constant "£eff~0.12 
below ~ 10 keV^". This statement is not correct. As 
explained and shown in Fig. 1 in the preprint, and also 
shown in Fig. 1 here, the £eff used in our analysis is not 
based on a single measurement (given the apparent dis- 
agreement within the quoted errors) but is the result 
from a global fit to all existing direct neutron scattering 
measurements. These results, where the recoil energy is 
measured by the scattering angle of mono-energetic neu- 
trons, are systematically cleaner than results that are 
based solely on a comparison between neutron calibra- 
tion data (featuring single nuclear recoils from a contin- 
uous neutron source) and Monte Carlo simulations. This 
is the reason why we have chosen not to include these 
£eff results in the global fit, and we have clarified this 
in the preprint. For illustration, we have added them in 
Fig. 1 of this note. There are two studies of this type 
(while only one is mentioned in [2]). The first is from 



the XENON 10 collaboration [5] and lies above our global 
fit. It is consistent with the 90% upper contour of £eff 
and clearly supports a constant extrapolation of >Ceff to 
lower energies. The ZEPLINIH results [6] do not agree 
with any of the tagged- neutron measurements, but we 
note that the assumed field quenching at the high oper- 
ating field of that experiment (3.9 kV/cm) has not been 
determined experimentally. 

Collar & McKinsey also point to a recent workshop 
presentation by a former Ph.D. student of the XENONIO 
collaboration (their ref. [8]) as "the most recent" analysis 
by the XENONIO collaboration. This result has not even 
been reviewed within the XENONIO collaboration and is 
unpublished. If this adds anything to this discussion, it 
underscores the point that the comparison of simulations 
and wide-spectrum neutron calibration data is subject to 
large systematic uncertainties. Furthermore, the lowest 
data points from the referenced slides were omitted in 
[2]. These points show a trend which is less supportive 
of their argument. 

We thus hold that at this stage an unbiased treatment 
of all the direct scattering data sets is the best approach 
and perform a global fit over all these data for the anal- 
ysis of XENONIOO data.This global fit >Ceff is monoton- 
ically falling with energy and its systematic uncertainty 
is given by the ±90% confidence level contours. By us- 
ing the correct statistical approach of a global fit with 
corresponding uncertainty, we have chosen not to favor 
our own £eff measurements over others. On the other 
hand. Collar & McKinsey suggest the use of one of the 
authors' own >Ceff measurement [7]. We do not agree with 
their statement that this "is featuring the best control of 
systematics so far". The trigger efficiency at the lowest 
energies probed in [7] is very low, and is not measured 
but inferred from a Monte Carlo simulation. The way 
this rapidly falling efficiency is corrected impacts sub- 
stantially the derived >Ceff- In the measurement of Aprile 
et al. [8], the trigger efficiency was >90% at 5 keV^, and 
this was directly measured and verified with a simula- 
tion. This means that all data points were obtained with 
similar systematics in [8] , a statement that does not hold 
for [7]. The systematic uncertainty from the spread in 
the scatter angle due to the finite size of the detectors is 
also larger in [7]. 

2. In their footnote, the authors express their confu- 
sion about the >Ceff curves used in the analysis and accuse 
us of trying to mislead the reader. We reject this notion. 
To avoid any misunderstandings, we have changed the 
text slightly to be more clear: In Fig. 1 of [1] (also in 
Fig. 1 here) there are 3 >Ceff contours: The thicker one in 
the center is the result from the global fit with a constant 
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FIG. 1: All published data on £eff- The black datapoints 
- used for the global fit in the XENONIOO paper [1] - are 
all published direct measurements of £eff- The red data 
(Sorensen (XENONIO) [5] and Lebedenko (ZEPLINIII) [6]) 
are from comparisons of data with Monte Carlo simulations. 
They were not used on the global fit because of their possibly 
larger systematic uncertainties. The three blue solid contours 
are the result from a global fit to all direct measurements 
(black) in the region from 5-100 keVr- The thinner contours 
above and below are the ±90% confidence level contours. The 
dashed lines below 5 keVr are the extrapolations as explained 
in the text. For the first XENONIOO data analysis, only the 
best fit and the lower 90% CL contour are used. 



extrapolation below 5 keV^ as explained above. The thin- 
ner contours above and below are the ±90% confidence 
level contours from this fit. To be very conservative, the 
lower contour is logarithmically extrapolated to energies 
below 5 keVr, with >Ceff = around 1 keV^. The slope 
of the extrapolation is far from "arbitrary" but fixed by 
a fit to the low energy part of the Yale points [7] and 
matched to the lower 90% confidence contour at 5 keV^. 

The logarithmic extrapolation is very conservative 
since a linear extrapolation describes the low energy part 
of the data points from ref. [7] equally well, and would re- 
sult in a much higher £eff and hence stronger constraints 
on low-mass WIMPs. From the three contours in Fig. 1, 
only the central ("global fit") and the lower one ("lower 
90% CL contour") are used in the XENONIOO analysis, 
as clearly stated in [1]. 

3. No satisfactory theory describing the behaviour of 
>Ceff in liquid xenon exists so far. The authors state that 
a kinematic cutoff to the production of scintillation is 
expected whenever the minimum excitation energy Eg of 
the system exceeds the maximum possible energy transfer 
to an electron by a slow-moving recoil ion, ^max- They 
refer to papers by Ahlen&Tarle [9] and Ficenec, Ahlen, 
Tarle et al. [10]. These papers deal with protons in or- 
ganic scintillators. Their arguments do not necessarily 
apply to Xe-Xe collisions. It is known in fact that Lind- 
hard theory [11, 12] is not adequate at very low energies, 
where mostly the tails of the ion-ion potential are probed 



and the Thomas-Fermi treatment becomes a crude ap- 
proximation. For Xe-Xe collisions this corresponds to 
about 10 keVr- The electron cannot be treated separately 
from the Xe atom and the maximum energy transferred 
to an electron cannot be given by simple kinematics, as 
advocated in [2]. 

The collision mechanism for heavy ions at very low 
energies may be better described by, e.g., the molecu- 
lar orbit theory [13], which involves many-body kinemat- 
ics. The argument by Collar and McKinsey is based on 
two-body kinematics and would not apply for heavy ion 
collisions in the energy region concerned here. In fact, 
Ficenec et al. [10] state that "No evidence for a response 
cutoff is observed at velocities extending well below the 
electron-excitation threshold of 6 x 10~^c expected from 
two-body kinematics" even for protons. Besides, if ^max 
for Xe-Xe is 39 keVr, the kinematics argument cannot 
explain the scintillation observed below 39 keV^ at all. 
Apart from the uncertainty in stopping power calcula- 
tions which affect directly nuclear quenching, other fac- 
tors may affect >Ceff through electronic quenching. How- 
ever, the current experimental and theoretical situation 
is such that there is no proven mechanism which justi- 
fies a decreasing i2eff with decreasing energy, as strongly 
advocated by Collar and McKinsey. 

We are fully aware of the impact of >Ceff on the over- 
all sensitivity of noble liquid dark matter experiments 
and our answer is simply that we will measure it again, 
extending it to the lowest possible energies. We need 
accurate data on this quantity and, within the XENON 
collaboration, we have already developed two new and in- 
dependent set-ups optimized to measure the energy and 
field dependence of both electron and nuclear recoils in 
liquid xenon. 

4. Finally, Collar and McKinsey doubt that we have 
properly taken into account the effects of the low number 
of photoelectrons at our threshold. (Note that this effect 
had not been accounted for in the preliminary plots pre- 
sented in their reference [17].) We agree that this has a 
crucial impact on the XENONIOO sensitivity to low mass 
WIMPs, however, it is a fact that an imperfect thresh- 
old due to a finite energy resolution leads to a mixing of 
events below threshold into the sample and vice versa. 
Since the expected WIMP spectrum is a steeply falling 
exponential (see Fig. 2), many more sub-threshold events 
fall in the energy region above threshold than vice versa. 

Due to the low number of detected photoelectrons at 
the XENONIOO threshold, the energy resolution is com- 
pletely dominated by counting statistics, therefore the 
expected true differential rate is convoluted with a Pois- 
son function to account for this behavior. We also point 
out that the XENONIOO efficiency is still very high down 
to 3 PE. 

Figure 2 shows the effect of Poisson broadening of our 
threshold for a DAM A benchmark case: There is a small 
amount of rate from a 10 GeV/c^ WIMP with a cross 
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FIG. 2: Expected spectrum of a 10 GeV/c^ WIMP with a 
cross section of 1 x 10~^^ cm^ (black, solid), a benchmark 
case at the lower edge of the DAMA region. The red (dashed) 
lines show the spectrum after a convolution with a Poisson 
distribution, the blue (thick dashed) line is corrected for the 
XENONIOO efficiency. The straight fines are the 3 PE and 
4 PE thresholds using the lower 90% CL £eff contour of the 
global fit as explained in the text. 
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FIG. 3: 90% confidence limits for the global fit of £eff with a 
threshold of 4 PE (black), and curves for the 90% lower con- 
tour of £eff at thresholds of 4 PE (yellow) and 3 PE (magenta). 
We used the following astrophysical parameters: galactic es- 
cape velocity = 544 km/s, WIMP density = 0.3 GeV/cm^, 
solar velocity = 220 km/s. 



section of 1 x IQ-"^^ cm^ leaking into the XENONIOO 
signal region, even at a threshold of 4 PE, correspond- 
ing to 9.6 keVr in the case of the lower 90% CL i2eff 
contour. Based on the light WIMP interpretation of the 
DAMA annual modulation signal, we would expect to see 
a total of 4.4 events above a lowered threshold of 3 PE 
(8.2 keVr for the conservative >Ceff case), taking into ac- 
count our reduced detection efficiency at lower energies 
(more than 18 events are expected for the best fit >Ceff 
above 3 PE corresponding to 7.0 keV^). As explained in 
our manuscript (see Fig. 3 in [1]) no event is observed 
leading to an exclusion of this case at 90% CL. 

The same effect at the lower end of the WIMP mass 
range favored by CoGeNT (7 GeV/c^ WIMP, cross sec- 
tions of 0.5 - 1 X 10~^^ cm^) reduces the expectation to 
0.73 - 1.5 events above a threshold of 3 PE in the case 
of the lower 90% CL £eff contour. Hence no significant 
conflict is found under this assumption and with the lim- 
ited exposure used so far. For the global flt with constant 
extrapolation of i2eff, this region is excluded at 90% CL. 

Fig. 3 shows the impact of the conservative choice of 
the 90% lower contour of £eff on the upper limits de- 
rived from the flrst XENONIOO dataset, both for our 
pre-deflned threshold of 4 PE and a lower threshold of 
3 PE. In general, other uncertainties play a role at low 
WIMP masses as well, e.g., the galactic escape velocity. 
We used a value of 544 km/s [14] . We further used a lo- 
cal WIMP density of 0.3 GeV/cm^ and a solar velocity of 
220 km/s, different (and less optimistic) from the values 
used in [2]. Our limits are distinctly different from the 
curves presented in [2]. 

In conclusion, we agree with the authors that the cur- 



rent situation on £eff in LXe is far from optimal and 
must be clarifled especially at the lowest Xe recoil ener- 
gies. However, in our manuscript we have properly taken 
into account the uncertainty by using >Ceff obtained from 
a global flt to all published direct measurements and by 
cross-checking the results with the lower 90% CL con- 
tour together with a very conservative extrapolation for 
Er <5 keVr- 
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